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ABSTRACT 

We present an updated spectroscopic orbit and a new visual orbit for the 
double-lined spectroscopic binary cr^ Coronae Borealis based on radial velocity 
measurements at the Oak Ridge Observatory in Harvard, Massachusetts and 
interferometric visibility measurements at the CHARA Array on Mount Wilson. 
(T^ CrB is composed of two Sun-like stars of roughly equal mass in a circularized 
orbit with a period of 1.14 days. The long baselines of the CHARA Array have 
allowed us to resolve the visual orbit for this pair, the shortest period binary yet 
resolved interferometrically, enabling us to determine component masses of 1.137 
± 0.037 Mq and 1.090 ± 0.036 Mq. We have also estimated absolute V"-band 
magnitudes of Mv(primary) = 4.35 ± 0.02 and Mv(secondary) = 4.74 ± 0.02. A 
comparison with stellar evolution models indicates a relatively young age of 1-3 
Gyr, consistent with the high Li abundance measured previously. This pair is the 
central component of a quintuple system, along with another similar-mass star, 

CrB, in a ~ 730-year visual orbit, and a distant M-dwarf binary, cr CrB C, 
at a projected separation of ~ 10'. We also present differential proper motion 
evidence to show that components C & D (ADS 9979C & D) listed for this 
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system in the Washington Double Star Catalog are optical alignments that are 
not gravitationally bound to the o CrB system. 

Subject headings: binaries: spectroscopic - stars: fundamental parameters - stars: 
individual (cr^ Coronae Borealis) - techniques: interferometric 



Introduction 



a CrB is a hierarchical multiple system 22 pc away. Its primary components, CrB (HR 
6064; HD 146362) and CrB (HR 606 3; HD 146361), are in a visual orbit with a preliminary 
period of ~ 900 years (jScardial Il979l ). of which the latt er is an RS CVn binary with a 



circularized and synchronized orbit of 1.139-day period (IStrassmeier fc Ricd l2003l . SR03 



hereafter). In addition to these three solar- type stars, the Washington Double Star Catalog^ 
(WDS) lists three additional com ponents for th is system. WDS components C and D wer e 
resolved 18" away at 103° in 1984 flPopovidl 19861 ) and 88" away at 82° in 1996 JCourtotlll996[ ). 
respectively. We will show in §[7] that both these components are optical alignments that are 
not gravitationally bound to the a CrB system. Finally, WDS com ponent E (a CrB C, HIP 
79551) which was resolved 635" away at 241° in 199 1 by Hipparco s ( jPerryman fc ESAlll997l ). 
was identified as a photocen tric-motion binary b v iHeintzl (jl990l ). The parallax and proper 
motion listed for this star in Ivan LeeuwenI (120071 ). the improved Hipparcos results based on 
a new reduction of the raw data, match the corresponding measures for cr^ CrB within the 
errors, confirming a physical association. 

SR03 presented photometric evidence in support of a rotation period of 1.157 ± 0.002 
days for both components of cx^ CrB, the central pair of this system. They explained the 
0.017- d ay dif ference between the rotation and orbital periods as differential surface rotation. 
Bakod (jl984j ) estimated an orbital inclination of 28°, assuming component masses of 1.2 
Mq based on spectral types. SR03 subsequently adopted this inclination to obtain com- 
ponent masses of 1.108 ± 0.004 Mq and 1.080 ± 0.004 Mq, but these masses are based 
on circular reasoning, and the errors are underestimated as they ignore the uncertainty in 



inclination. Several spectroscopic orbit s have been published for this pair (IHarperl Il925 



Bakoslll984l : iDuquennoy fc Mayorlll99ll : SR03), enabling the spectroscopic orbital elements 
to be well-constrained. We present an updated spectroscopic solution based on these prior 
data and our own radial velocity measurements (§ 12. 1^ § 14. ip . Our visual orbit leverages 
these spectroscopic solutions and derives all orbital elements for this binary (§E]), leading to 



http : / /ad . usno . navy, mil / wds / 
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accurate component masses (§ l6.ip . 



This work uti l izes a very precise parallax measure for this radio-emitting binary obtained 
by lLestrade et al.l (119991 ) using Very Long Baseline Interferometry (VLBI). Their parallax of 
43.93 ± 0.10 mas is about 10 times more preci se than the Hipparc os catalog value of 46.11 
±0.98 m as and 12 times more precise than the Ivan LeeuwenI (120071 ) measure o f 47.35 ± 1.20 
mas. The iLestrade et al.l value is 2.2-cr and 2.9-0" lower than the Hipparcos and Ivan Leeuwen 
measures, respectively. To check for systematic offsets, we compared the parallaxes for all 
overlapping stars in these three sources. While the differ ence in parallax is most significant for 
(T^ CrB, we found no systematic differences. Moreover, iLestrade et al.l performed statistical 
checks to verify the accuracy of their measure, so we adopt their parallax to derive the 
physical parameters of the component stars (§[!])• 

The Center for High Angular Resolution Astronomy (CHARA) Array's unique capa- 
bilities, facilitated by the world's l ongest optical interferometric baselines, have enabled a 



varie ty of astrophysical studies (e.g.. lMcAlister et al.ll2005l : iBaines et al.ll2007l : iMonnier et al. 



20071 ). This work utilizes the Array's longest baselines to resolve the 1.14-day spectroscopic 
binary, the shortest period system yet resolved. While this is the first visual orbit determined 
using interferometric visibilities measured with the CHARA Array, the technique described 
here has regular ly been employed fo r longer-period bin aries using other long-baseline inter- 



ferometers (e.g.. lHummel et al.lll993l : lBoden et al.lll999l ). The CrB binary has a projected 



angular separation of about 1.1 mas in the sky, making it easily resolvable for the CHARA 
Array, which has angular resolution capabilities in the K' band down to about 0.4 mas for 
binaries. 



2. Spectroscopic Measurements 

Spectroscopic observations of cr^ CrB were conducted at the Harvard-Smithsonian Cen- 
ter for Astrophysics (CfA) with an echelle spectrograph on the 1.5m Wyeth reflector at the 
Oak Ridge Observatory in the town of Harvard, Massachusetts. A total of 46 usable spectra 
were gathered from 1992 May to 1999 July, each of which covers a single echelle order (45 A) 
cent ered at 5188.5 A and was recorded using an intensified photon-counting Reticon detector 



see 



Lathamlll992l ). The strongest lines in this window are those of the Mg I b triplet. The 
resolving power of these observations is A/AA ~ 35,000, and the nominal signal-to-noise 
ratios range from 21 to 94 per resolution element of 8.5 kms~^. 



Radial velocities were obtain ed using the two-dimensional cross-correlation algorithm 
TODCOR (IZucker fc Mazehlll994j ). Templates for the cross correlations were selected from 
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an e xtensive library of calculated spectra Tjased on model atmospheres by R. L. Kuruc^ (see 
also iNordstrom et al.lll994j : iLatham et al.ll2002l ). These calculated spectra cover a wide range 
of effective temperatures (Teg), rotational velocities [vsini when seen in projection), surface 
gravities (log (7), and metallicities. Experience has shown that radial velocities are largely 
insensitive to the surface gravity and metallicity adopted for the templates. Consequently, 
the optimum template for each star was determined from extensive grids of cross-correlations 
varying the temperature and rotational velocity, seeking to maximize the average correlation 
weighted by the strength of each exposure. The results we obtain, adopting logg = 4.5 and 
solar metallicit}|f] for both stars, are T^g = 6050 K and vsini = 26 kms"^ for the primary, 
and Tefj = 5870 K and vsini = 26 kms~^ for the secondary. Estimated uncertainties are 
150 K and 1 kms~^ for the temperatures and projected rotational velocities, respectively. 
Template parameters near these values were selected for deriving the radial velocities. Typ- 
ical uncertainties for the velocities are 1 kms~^ for both stars. 

The stability of the zero-point of our velocity system was monitored by means of ex- 
posures of the dusk and dawn sky , and small run-to-run corrections were applied in the 



19921). Additi o nal co rrecti ons for systematics were applied to 



Latham et al.1 (119961 ) and [Torres et al.l (119971 ) to account for 



manner described by [Latham 
the velocities as described by 
residual blending effects. These corrections are based on simulations with artificial composite 
spectra processed with TODCOR in the same way as the real spectra. The final heliocentric 
velocities and their l-cx errors are listed in Table [H along with the corresponding epochs of 
observation, — C residuals, and orbital phase. 

The light ratio between the components was estimated directly from the spectra follow- 
ing |Zu£kei_&.^:Iazeh| (119941 ). After corrections for systematics analogous to those described 
above, we obtain = 0.67±0.02 at the mean wavelength of our observations (5188.5 A). 
Given that the stars have slightly different temperatures, a small correction to the visual 
band was determined from synthetic spectra integrated over the V passband and the spectral 
window of our observations. The corrected value is (4/^?);/ = 0.70 ± 0.02. 



The visual companion CrB was also observed spectroscopically at the CfA with the 



^Available at http: //cf aku5 . cf a. harvard. edu| 

•^SROS ha ve reported a me t allicit y for cr^ CrB of [Fe/H] = —0.37 with an uncertainty no smaller than 
0.1 dex, and INordstrom et all (|2004f ) reported the value [Fe/H] = —0.24 based on Stromgrcn photometry. 
Metallicity determinations for double- lined spectroscopic binaries are particularly difficult, and both of these 
estimates are likely to be affected at some level by the double-lined nature of the system. However, the visual 
co mpanion (a^ CrB) i s appa rently a single star, and has an accurate spectro scopic abundance determination 
bv lValenti fc Fischer] |2005h giving [Fe/H] = -0.06 ± 0.03, and another bv lFuhrmaniil |2004l) giving [Fc/H] 
= — 0.064±0.068. The near-solar metallicity from these determinations is considered here to be more reliable. 
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same instrumental setup. We obtained 18 observations between 1996 June and 2004 August. 
The stellar parameters were determined with a procedure similar to that used for cr^ CrB, 
and yielded Tgfj = 5950 ± 100 K and vsini = 3 ± 2 kms~^, for an adopted logg = 4.5 
and solar metallicity (see Footnote [3]). Radial velocities were obtained with standard cross- 
correlation techniques using a template selected according to the above parameters. These 
measurements give an average velocity of — 14.70±0.11 kms~^, with no significant variation 
within the observational errors. We use this radial velocity to unambiguously determine the 
longitude of the ascending node for the wider a^ — a"^ CrB visual orbit (§ 16.41) . 



2.1. Historical Data Sets 



In additio n to our own , four other r adial- velocity data sets hay e been published in 



the literature (IHarperl Il925l : iBakod Il984l : iDuquennoy fc Mayorl Il99ll : SR03). Except for 



the more recent one, the older data are generally of lower quality and contribute little to 
the mass determinations, but they do extend the time coverage considerably (to nearly 86 
years, or 27,500 orbital cycles) and can be used to improve the orbital period. Because of 
our concerns over possible systematic differences among different data sets, particularly in 
the velocity semi-amplitudes but also in the velocity zero points, we did not simply merge 
all these observations together indiscriminately, but instead we proceeded as follows. We 
considered all observations simultaneously in a single least-squares orbital fit, imposing a 
common period and epoch of maximum primary velocity in a circular orbit, but we allowed 
each data set to have its own velocity semi- amplitudes {Kp, Kg) as well as its own systematic 
velocity zero-point offset relative to the reference frame defined by the CfA observations. 
Additionally, we included one more adjustable parameter per set to account for possible 
systematic differences between the primary and secondary velocities in each group. These 
were statistically significant only in the observations by SR03. Relative weights for each data 
set were determined by iterations from the RMS residual of the fit, separately for the primary 
and secondary velocities. The resulting orbital period is P = 1.139791423 ± 0.000000080 
days, and the time of maximum primary velocity nearest to the average date of the CfA 
observations is T = 2,450,127.61845 ± 0.00020 (HJD). We adopt this ephemeris for the 
remainder of the paper. 



3. Interferometric Measurements 



Interferometric visibilities for o"^ CrB were measured during 2007 May— July at the 
CHARA Array's six-element long-baseline interferometer located on Mount Wilson, Call- 
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fornia (jten Brummelaar et al.ll2005l ). The Array uses the visible wavelengths 480-800 nm 
for tracking and tip/tilt corrections, and the near-infrared K' (2.13 /im) and H (1.67 yum) 
bands for fringe detection. The 26 visibility measurements used in the final orbit determina- 
tion, listed in Table [21 were obtained in the K' band on the Sl-El and S1-E2 two-telescope 
baselines spanning projected baselines of 268-331 meters. The interference fringes were ob- 
tained using the pupil-plane "CHARA Classic" beam combiner. While some of the data were 
obtained via on-site observing at Mount Wils on, the bulk of the data were gathered at the Ar- 
rington Remote Operations Center (AROC. iFallon et al.ll2003l ) located at the Georgia State 
University campus in Atlanta, Georgia. Following the standard practice of time-bracketed 
observations, we interleaved each target visibility measurement with those of a calibrator star 
(HD 152598) in order to remove instrumental and atmospheric e ffects. For furth e r deta ils 
on the observing practice and the data reduction process, refer to iMcAlister et al.l (120051 ). 

We selected HR 6279 (HD 152598), an FOV star offset from CrB by 8?3, as the 
calibrator based on its small estimated angular diameter and its apparent lack of any close 
companions . We obtained phot omet ric measurements for th is star in the Johnson UBV 
bands from iGrenier et al.l (1 19851 ) and IPerryman &: ES7\I (119971 ) , and JHKs bands from the 
Two Micron All Sky Surveju (2MAS S) and transformed ther n to calibrat e d flux measure- 
ments using the methods described in lColina et al.l (119961 ) and lCohen et al.l (120031 ). We then 
fitted these fluxes to spectral energy distribution model^, yielding an angular diameter of 
0.467 ± 0.013 mas for HD 152598, corresponding to T^e = 7150 K and logg = 4.3. This 
diameter estimate results in a predicted calibrator visibility of Kai = 0.858 ± 0.008 at our 
longest baseline of 330 m, contributing roughly 1% error to the calibrated visibilities. This 
error is included in our roughly 10% total visibility errors listed in Table [21 along with the 
epoch of observation (at mid-exposure), the target star's calibrated visibility, the predicted 
visibility for the best-fit orbit, the — C visibility residual, the baseline projections along 
East- West (u) and North-South {v) directions, and the hour angle of the target. 



4. Determination of the Orbit 

Consistent with prior evidence of a synchronized orbit (SR03), we adopt a circular orbit 
(e = 0, u; = 0) with the orbital period (P) and epoch of nodal passage (T) from § l2.1l for the 
spectroscopic and visual orbit solutions presented below. 



^http : //www . ipac . caltech . ed.u/2inass 

^The model fliaxes were interpolated from the grid of models from R. L. Kurucz, available at 
http : //cf aku5 . cf a . harvaxd . edu 
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4.1. Spectroscopic Orbital Solutions 

Our measured radial velocities enable us to derive the three remaining spectroscopic or- 
bital elements, namely, the center-of-mass velocity (7) and the radial velocity semi-amplitudes 
of the primary and secondary {Kp and Kg, respectively). To check for consistency with prior 
efforts, we used the velocities published in SR03 to derive a second orbital solution. The 
calculated radial velocities for the derived orbits are shown in Figures [1] and [2] (solid and 
dashed curves for the primary and secondary, respectively) along with the measured radial 
velocities and residuals for the primary (filled circles) and secondary (open circles). The 
corresponding orbital solutions are presented in Table [3] along with the related derived quan- 
tities. For comparison purposes, we have also included the values presented in SR03, which 
are consistent with our orbit generated using their velocities. However, the orbit obtained 
using our velocities is statistically different from the one obtained using SR03 velocities. 
While the primary's velocity semi-amplitude matches within the errors between these two 
solutions, the secondary's differs by over 5-cr, resulting in a 4-cr difference in the mass ratios. 

One possible explanation of the difference in the orbital solutions could be the velocity 
residuals for the orbit using SR03 data (Figure [2]), which show an obvious pattern for both 
components. Those observations were obtained on four nights over a five-day period. To 
further examine these patterns, we display the residuals for each of the four nights in Figure [31 
as a function of time. Clear trends are seen on each night, which are different for the primary 
and secondary components and have peak-to-peak excursions reaching 4 kms~^ in some 
cases, significantly larger than the velocity errors of 0.1-1.2 kms~^ (SR03). On some but 
not all nights, there appears to be a periodicity of roughly 0.20-0.25 days. The nature of 
these trends is unclear, particularly because this periodicity is much shorter than either the 
orbital or the rotational periods. Instrumental effects seem unlikely, but an explanation in 
terms of the considerable spottedness of both stars is certainly a distinct possibility. The 
Doppler imaging maps produced by SR03 show that both components display a very patchy 
distribution of surface features covering the polar regions. Individual features coming in and 
out of view as the stars rotate could easily be the cause of the systematic effects observed 
in the radial velocities, and the effects would not necessarily have to be the same on both 
stars, just as observed. Slight changes in the spots from one night to the next could account 
for the different patterns seen in Figure [31 The relatively large amplitude of the residual 
variations raises the concern that they may be affecting the velocity semi-amplitudes of the 
orbit, depending on the phase at which they occur. We do not see such trends in the CfA 
data, perhaps because our observations span a much longer time (more than 7 years, and 
^^2200 rotational cycles), allowing for spots to change and average out these effects. We 
therefore proceed on the assumption that possible systematic effects of this nature on Kp 
and Ks are lessened in the CfA data. 
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5. The Visual Orbit Solution 

The basic measured quantity from an interferometric observation is visibility, which 
evaluates the contrast in the fringe pattern obtained by combining starlight wave fronts 
from multiple apertures, filtered through a finite bandwidth. For a single star of angular 
diameter 6, the interferometric visibility V for a uniform disk model is given by, 

_ 2J,inBe/X) 

ttBO/X ' ^ ' 

where Ji is the first-order Bessel function, B is the projected baseline length as seen by the 
star, and A is the observed bandpass central wavelength. The interferometric visibility for 
a binary, where the individual stars have visibilities Vp (primary) and Vs (secondary) per 
equation ([1]) above, is given by 

JW%' + K' + 2/3V;V.cos((2t/A)B . s)) 
—3 • 

where f3 is the primary to secondary flux ratio, B is the projected baseline vector as seen by 
the binary, and s is the binary's angular separation vector in the plane of the sky. 

Using our measured interferometric visibilities and the above equations, we are able to 
augment the spectroscopic orbital solutions to derive a visual orbit for cr^ CrB. Adopting the 
period and epoch of nodal passage from § 12. H we now derive the parameters that can only 
be determined astrometrically: angular semimajor axis (a); inclination (i); and, longitude of 
the ascending node {Q). We also treat the i^'-band magnitude difference as a free parameter 
in order to test evolutionary models. 

For a circular orbit, the epoch of periastron passage (Tq) is replaced by the epoch of 
ascending nodal passag e (T„ade), def ined as the epoch of fastest secondary recession, in the 



visual orbit equations (IHeintzl Il978l ). Accordingly, we translate the T value listed in § 12.11 
by one-half of the orbital period to determine the epoch of the ascending nodal passage as 
Tnode = 2,450,127.04855 ±0.00020 (HJD) for use in our visual orbit solution. The 1-a errors 
of this and other adopted parameters listed in Table H] have been propagated to our error 
estimates for the derived parameters. 

The angular diameters of the components are too small to be resolved by our i^^'-band 
observations. We therefore estimate these based on the components' absolute magnitudes 
and temperatures as described below. We first estimate the Johnson V^-band magnitude of 
CrB using its Tycho-2 magnitudes of B^ = 6.262 ± 0.014 and Vt = 5.620 ± 0.009 and 
the relation Vj = Vt — 0.090(i?T — Vr) from the Guide to the Tycho-2 Catalog. Then, using 
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the \^-band flux ratio from §|2] and tlie iLestrade et al.l (Il999l ) parallax, we obtain absolute 
magnitudes of My = 4.35 ± 0.02 for the primary and My = 4.74 ± 0.02 for the secondary. 
These magnitudes lead to linear radius estimates of 1.2i?0 for the prin aary and l.l-R^ for 



the secori d ary u sing the tabulation of stellar physi c al pa rameters in iPopperl (jl980l ) and 
AndersenI ( 1l99ll ). Finally, using the ILestrade et al.l (Il999l ) parallax, we adopt component 
angular diameters of 9p = 0.50 mas and = 0.45 mas, propagating a 0.05 mas uncertainty 
in these values for deriving the uncertainty of our orbital elements. Diameter estimates using 
the temperatures of the components from §|2] are consistent with these values. 

We conduct an exhaustive search of the parameter space for the unknown parameters 
mentioned, namely, a, i, Q, and AK'. The orb ital inclination i s con strained by the asini 
from spectroscopy, the free-parameter a, and the ILestrade et al.l (119991 ) parallax. We impose 
this constraint during our exploration of the parameter space along with its associated 1-a 
error. We explore the unknown parameters over many iterations, by randomly selecting them 
between broad limits and using equation ([2]) to evaluate the predicted binary visibility for 
the baseline and binary positions at each observational epoch. The orbital solution presented 
here represents the parameter set with the minimum value when comparing the predicted 
and measured visibilities. 

Figure H] shows the measured visibilities (plus signs) with vertical error bars for each 
of the 26 observations, along with the computed model visibilities (diamonds), and Table [2] 
lists the corresponding numerical values of the observed and model visibilities along with 
the residuals of the fit. Table H] summarizes the visual orbit parameters for cr^ CrB from our 
solution and Figure [5] plots the visual orbit in the plane of the sky. As seen in Figure |5l we 
have a reasonably good phase coverage from our observations. 



As mentioned in § 14.11 star spots can create systematic effects in the data obtained 
on this binary. These effects are especially significant for data obtained over a short time 
baseline, as seen for the SR03 spectroscopic solution. While our interferometric data span 
73 days, allowing for some averaging of these effects, the bulk of the data used were obtained 
over 12 days, justifying an exploration of this effect. Specifically, the separation between the 
stars derived from our visibility data woul d represent the separ ation of the centers of light 
rather than that of mass. As discussed in iHummel et al.l (119941 ). heavily-spotted stars will 
incur a systematic shift in the center of light from rotational and orbital motions, perhaps 
inducing an additional uncertainty in the orbital elements derived. We assume a spot-induced 
change in the angular semimajor axis of 2% of the primary's diameter, or 0.01 mas. This 
is less that the uncertainty of our derived semimajor axis, and at our baselines of 270-330 
meters, translates to a 0.005-0.011 change in the visibility. While the uncertainties of our 
measured visibilities are an order of magnitude larger than this, we ran a test orbital fit by 
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adding a 0.010 uncertainty to the visibility errors as a root-sum-squared. While, as expected, 
the of the fit improved, the values and uncertainties of the derived parameters remained 
unchanged, leading us to conclude that this effect, while real, is too small to affect our 
results. 

We determine the 1, 2, and 3-cr uncertainties of each visual orbit parameter using a 
Monte Carlo simulation approach. We compute the orbital fit for 100,000 iterations, where 
for each iteration, we randomly select the adopted parameters within their respective 1-a 
intervals and the model parameters around their corresponding best-fit solution, generating 
a multi-dimensional "surface" . Then, we project this surface along each parameter axis, 
resulting in the plots shown in Figures [6] to [9l The figures show the distribution around the 
best-fit orbit and enable estimation of 1, 2, and 3-cr errors for each parameter based on a 
deviation of 1, 4, and 9 units, respectively, from its minimum value. The horizontal dashed 
lines in the figures from bottom to top mark the minimum value and those corresponding 
to 1, 2, and 3-cr errors, and Table Hllists the corresponding numerical 1-cr errors of the model 
parameters. 



6. Physical Parameters 
6.1. Component Mass Estimates 

Our angular semimajor axis obtained from i nterfe rometry translates to 0.0279 ± 0.0003 



AU or 5.99 ± 0.07 Rq using the iLestrade et al.l (Il999l ) parallax. Newton's generalization of 



Kepler's Third Law then yields a mass-sum of 2.227 ± 0.073 Mq for the pair, and using the 
mass ratio from our spectroscopic solution of 0.9586 ± 0.0047, we get individual component 
masses of 1.137 ± 0.037 Mq and 1.090 ± 0.036 Mq for the primary and secondary, respec- 
tively. As noted in § 14. H the SR03 velocities yield a significantly different mass ratio of 0.9746 
± 0.0016, but this 4-cr difference is not enough to influence the mass estimates significantly. 
The uncertainty in our masses is dominated by the cubed semimajor axis factor in estimating 
the mass sum, resulting in about a 3% uncertainty in mass-sum corresponding to a 1% un- 
certainty in the semimajor axis. The high precision of the mass ratio from the spectroscopic 
solution results in final masses of 3% uncertainty as well. Component mass estimates using 
the SR03 velocities are 1.128 ± 0.037 and 1.099 ± 0.036, in excellent agreement with the 
masses using our velocities. These masses along with other physical parameters derived are 
listed in Table [51 
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6.2. Radii of the Components 



Assuming synchronous and co-aligned rotation of spherical components, reasonable 
given the short orbital period and evidence from SR03 of unevolved stars contained within 
their Roche limits, we can estimate the component radii from the measured spectroscopic 
vsini. As mentioned in §|21 our spectra yield vsini = 26 ± 1 kms^^ for both the primary 
and secondary. These values and uncertainties are identical to those in SR03. Using the 
inclination from our visual orbit, and adopting the orbital period from spectroscopy as the 
rotational period, we get identical component radii of 1.244 ± 0.050 Rq for the primary and 
secondar y. This translates to a n angular diameter of 0.509 ± 0.020 mas for each component 
using the iLestrade et al.l (119991 ) parallax, in excellent agreement with our adopted diameter 
for the primary and a 1-a variance for the secondary, given our associated 0.05-mas errors 
for these values. These radii estimates, along with the effective temperatures from §[2] and 
the relation L (x TJ g, lead t o a lum inosity ratio of 0.89 ± 0.16. Alternatively, using bolo- 



metric corrections from lFloweii (119961 ) of BCp = -0.038 ± 0.017 and BC^ = -0.064 ± 0.020 



corresponding to the components' effective temperatures, the ^-band flux ratio of 0.70 ± 
0.02 from spectroscopy translates to a total luminosity ratio of 0.68 ± 0.20, a 1-a variance 
from the estimate above. Conversely, our estimates of effective temperature and luminosity 
ratio require a radius ratio of 0.88 ± 0.14, again at a 1-a variance from the 1.00 ± 0.06 
estimate from the identical v sin i values of the components. 



6.3. Absolute Magnitudes and Ages 



We allowed the i^'-band magnitude difference to be a free parameter for our visual orbit 
fit, obtainin g AK' = 0.19 ±0.19 , consi stent with the 0.18 estimate from the mass-luminosity 
relations of 



Henrv fc McCarthvl Jl993l lR. The uncertainty in AK' is large because visibility 
measurements of nearly equal mass, and hence nearly eq ual brightness, pairs are relatively 



insensitive to the magnitude difference of the components (iHummel et al.lll998l : iBoden et al. 
19991 ). Using equation ([2]), we have verified that a 10% change in AK' for CrB results 
in only 0.1% change in visibility. This, along with the poor quality K magnitude listed in 
2MASS (for o"^ CrB, K = 4.052 ±0.036, but flagged as a very poor fit), thwart any attempts 
to use these magnitudes for checking stellar evolution models. However, we can revert to 
y-band photometry to explore this topic. 



^The relations from iHenrv fc McCarthvl are for 0.5 Mq < Mass < 1.0 Mq. We consider it safe to 
extrapolate out to our estimated masses of slightly larger than 1.0 Mq. 
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In §[5], we derived the absolute ^-band magnitudes of the components of CrB as 
My = 4.35 ± 0.02 for the primary and My = 4.7 4 ± 0.02 for th e secon dary. For CrB, 
we similarly use the Tycho-2 magnitudes and the iLestrade et al.l (119991 ) parallax to obtain 
My = 4.64 ± 0.01. SR03 had a smaller magnitude difference for the components of cr^ CrB, 
and the corresponding results using their spectroscopy are also included in Table [5] along with 
the values from their paper. Figure [10] plots these three stars on an H-R diagram using our 
magnitude and temperature estimates, along with isochro nes for 0.5, 1.5 , 3.0, and 5.0 Gyr 
ages (left to right) from the Yonsei-Yale is ochrones (dotted, lYi et al.ll2001l) and the Victor ia- 
Regina stellar evolution models (dashed, VandenBerg et al. 20061 ) for solar metallicity (see 
Footnote [3]). 



Wright et al.l (120041) estimate an age of 1.8 Gyr for CrB based on chromospheric 



activity, and IValenti fc Fischerl (120051 ) estimate an age of 5.0 Gyr from spectroscopy with 
limits of 2.9-7.8 Gyr based on l-cr changes to logL. SR03 identify a much lower age, of a 
few times 10^ years, by matching pre- main-sequence evolutionary tracks and point to their 
higher Li abundance as supporting evidence. While abundance determinations in double- 
lined spectroscopic binaries are particularly difficult and more prone to errors, the high Li 
abundance of 2.60 ± 0.03 (SR03) for the slow-rotating single-lined companion CrB does 
argue for a young system. Each point along the isochrones plotted in Figure [10] corresponds 
to a particular mass, allowing us to use our mass estimates for the components of CrB 
to further constrain the system's age. Our mass, luminosity, and temperature estimates 
indicate an age for this system of 0.5-1.5 Gyr, with a range of 0.1-3 Gyr permissible within 
1-0" errors. 



6.4. Mass Estimate of CrB 



Our mass estimates for the components of C r B allo w us to constrain the mass of 
the wider visual companion CrB as well. IScardial (1 1979! ) presented an improved visual 
orbit for the AB pair based on 886 observations spanning almost 200 years of observation, 
yielding P = 889 years, a = 5'.'9, i = 31?8, e = 0.76, and Q = 16?9. However, he did not 
publish uncertainties for these parameters, and given the long period, his less than one-third 
phase coverage leads to only a preliminary orbital solution, albeit one that convincingly 
shows orbital motion of the pair. He further uses parallaxes available to him to derive a 
mass-sum for the AB system of 3.2 MfT^. W e used all current WDS observations, adding 
almost 200 observations since IScardial (Il979l ). to update this orbit and obtain uncertainties 
for the parameters. Our visual orbit is presented in Figure [11], along with th e IScardial orbit 



for comparison, and Table [6] lists the derived orbital elements. Adopting the ILestrade et al. 
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( 119991 ) parallax of the A component, we estimate a mass-sum of 3.2 ± 0.9 Mq, re sulting in a 



B-component mass estini a te of 1.0 Mq, consistent with its spectral type of Gl V (iGray et al. 



20031 ) . IValenti fc Fischer! (120051 ) estimate a mass of 0.77 ± 0.21 Mq based on high-resolution 



spectroscopy, but we believe that they systematically underestimate their uncertainty by 
overlooking the log e factor in converting from uncertainty in log L to uncertainty in L. Using 
the loge factor, we followed their methods to obtain a mass estimate of 0.7 7 ± 0.44 Mq. 



Th e mass-error is dominated by the uncertainty of the lGliese fc Jahreifil (11991) para llax used 



by IValenti fc Fischen (120051 ) . Adopting the higher precision iLestrade et al. 



19991 ) parallax 



of the primary, we follow their method, and using the loge factor, get a mass estimate of 
0.78 ± 0.11 Mq. This mass is too low for the spectral type (as well as our own estimate 
of the effective temperature; see §[2]) and the expectation from the visual orbit. A possible 
contamination of the secondary's spectral type from the 7" distant primary is unlikely, as 
determined by Richard Gray at our request from new spectroscopic observations (R. Gray 
2008, private communication). 

The inclination and longitude of the ascending node for this visual orbit are similar to 
those of the inner (a^ CrB) orbit, suggesting coplanarity. For the outer visual orbit, we can 
use our radial velocity estimate for CrB, our derived systemic velocity for CrB, and 
the speckle observations to unambiguously determine the longitude of the ascending node 
as fl = 28° ± 0?5. Using the equation for the relative inclination of the two orbits (0) from 
Fekell ( 1l98ll ). we get = 4?7 or 60?3, given the 180° ambiguity in Q for the inner orbit. 



confirming coplanarity possibility. 



7. The Wide Components: Optical or Physical? 



In addition to the three solar-type stars, the WDS lists three additional components for 
a CrB. We present evidence to show that WDS components C and D are optical alignments, 
while component E, itself a binary, is a phy sical associatio n. WDS component C (ADS 
9979C), measured 18" away at 103° i n 1984 (Popovi^ 1986 ) has a proper motion of /ia = 



O'.'Oie yr-i and 
of cr^ CrB of Hn 



-0'.'015 yr-i Jjeffers et all Il963l ) 



van 



-0'.'26364 ± 0'.'00091 yr'^ and fis 



significantly different from that 
-0'.'09259 ± 0'.'00129 y r"^ from 



eeuweru (120071 ). Similarly, component D, measured 88" away at 82° in 1996 (jCourtot 



19961 ) and clearly seen by us as a field star by blinking the multi-epoch STScI Digitized Sky 
Surve:vlll (DSS) im ages, has a proper motion of /Iq, = +0'.'004 yr~^ and fxs = — 0'.'017 yr~^ 



( iJeffers et al.lll963l ). again significantly different from that of CrB. As a confirmation of the 



^ ,http : //stdatu . stsci . edu/cgi-bin/dss 



form 
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optical alignment, we compare in Figures [12] and [13] the observed separations of components 
C and D, respectively, from the primary with the corresponding expected values based on 
their proper motions. The solid line is a linear fit to the published measurements from the 
WDS and the dashed line is the expected separation based on differential proper motion. 
The excellent agreement between the two lines for both components confirms them as field 
stars. 

WDS component E (cr CrB C, HIP 79551) is widely separated from the p rimary at 635", 



transl ating to a minimum physical separation of over 14,000 AU using the iLestrade et al. 



( I1999I ) parallax. Despite its wide separation, this component appears to be physically as- 
sociated with a CrB based on its matching parallax of vr = 45.40 ± 3.71 m as and proper 
motio n of /i^ = -0'.'26592 ± 0'.'00299 yr'^ and fis = -0'.'08363 ± 0'.'00368 yr'^ dvan Leeuwen 
20071 ). While seemingly extreme for gravitationally bound systems, phy sical association has 



been demonstrated for pairs with separations out to 20,000 AU (e.g., [Latham et a. 



1991 



Poveda et al.[ [l994i ). a CrB C has a spectral classi fication of M2.5V (iReid et al.l [1995[ ). 
apparent magnitude of V = 12.24 (IBidelmanl [l985[ ). and has itself been identifie d as a 



photo c entric-motion bina r y wit h an unseen companion of 0.1 Mq in a 52-year orbit (iHeintz 
1990|). [Ferryman fc ESA[ (119971 ) also identifies this star as a binary of type 'X' or stochastic 
solution, implying a photocenter wobble for an unresolved star, but for which the Hipparcos 
data are not sufficient to derive an orbit. 



8. Conclusion 



Augmenting our radial velocity measurements with published values, we obtain a cov- 
erage of nearly 86 years or 27,500 orbital cycles, resulting in a very precise ephemeris of 
P = 1.139791423 ± 0.000000080 days and T = 2,450,127.61845 ± 0.00020 (HJD) and a ro- 
bust spectroscopic orbit for o"^ CrB. Using the CHARA Array, we have resolved this 1.14-day 
spectroscopic binary, the shortest period system yet resolved, and derived its visual orbit. 
The resulting component masses are 1.137 ± 0.037 Mq and 1.090 ± 0.036 Mq for the primary 
and secondary, respectively. Our spectroscopy supports prior efforts in estimating the same 
V sin i values for both components, which assuming a synchronized, co-aligned rotation, re- 
sults in equal radii of 1.244 ± 0.050 Rq for both components. The corresponding radius ratio 
is consistent within 1-a with its estimate using the components' temperatures and flux ratio 
from spectroscopy. We have also shown that this binary resides in a hierarchical quintuple 
system, composed of three close Sun-like stars and a wide M-dwarf binary. The wider visual 
orbit companion, CrB, is about 7" away in a 726-year visual orbit with i = 32?3, which 
appears to be coplanar with the inner orbit. A comparison of the mass and absolute magni- 
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tude estimates of CrB and CrB with current stellar evolution models indicates a young 
age for the system of 1-3 Gyr, consistent with the relatively high Li abundance previously 
measured. Finally, the widest member of this system is an M-dwarf binary, a CrB C, at a 
minimum separation of 14,000 AU. Figure [H] depicts the system's hierarchy in a pictorial 
form. 
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Table 1. Radial Velocities of ct^ CrB 



HJD RVp RV, cjRv^ <jRv^ (O-C)p {O - C), Orbital 
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Table 1 — Continued 
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Table 2. Interferometric Visibilities for cr^ CrB 
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(m) 


(h) 


54237.763 


0.864 


0.086 


0.783 


0, 


.081 


202, 


.4 


250.7 


-2.24 


54237.774 


0.909 


0.107 


0.775 


0, 


.134 


196, 


.7 


258.2 


-1.99 


54237.784 


0.736 


0.062 


0.759 


-0, 


.022 


190, 


.3 


265.2 


-1.74 


54237.796 


0.702 


0.063 


0.729 


-0, 


.027 


182, 


.4 


272.6 


-1.46 


54237.806 


0.585 


0.058 


0.688 


-0, 


.103 


174, 


.6 


278.9 


-1.22 


54237.816 


0.652 


0.076 


0.625 


0, 


.027 


165, 


.6 


285.3 


-0.97 


54237.833 


0.468 


0.053 


0.474 


-0, 


.006 


149, 


.7 


294.7 


-0.56 


54237.932 


0.833 


0.049 


0.833 


0, 


.001 


30, 


.4 


326.9 


1.82 


54237.942 


0.775 


0.059 


0.791 


-0, 


.017 


17, 


.1 


327.7 


2.05 


54237.954 


0.672 


0.038 


0.672 


0, 


.001 


0, 


.5 


328.1 


2.34 


54237.980 


0.244 


0.015 


0.247 


-0, 


.004 


-35, 


.3 


326.5 


2.98 


54247.701 


0.858 


0.113 


0.887 


-0, 


.029 


159, 


.9 


214.9 


-3.08 


54247.716 


0.888 


0.080 


0.863 


0, 


.025 


154, 


.1 


223.0 


-2.73 


54247.729 


0.824 


0.083 


0.785 


0, 


.040 


147, 


.5 


230.2 


-2.40 


54247.744 


0.669 


0.093 


0.644 


0, 


.025 


139 


,1 


237.6 


-2.05 


54247.761 


0.435 


0.058 


0.430 


0, 


.005 


128, 


.1 


245.6 


-1.64 


54249.714 


0.589 


0.053 


0.621 


-0, 


.032 


152, 


.1 


225.3 


-2.63 


54249.726 


0.570 


0.054 


0.609 


-0, 


.039 


146, 


.6 


231.1 


-2.36 


54249.739 


0.575 


0.064 


0.573 


0, 


.002 


138, 


.6 


238.1 


-2.03 


54249.751 


0.594 


0.063 


0.524 


0, 


.070 


131, 


.3 


243.5 


-1.75 


54249.772 


0.391 


0.059 


0.376 


0, 


.015 


115, 


.7 


252.8 


-1.24 


54310.716 


0.616 


0.062 


0.526 





.090 


48, 


,7 


325.0 


1.49 


54310.726 


0.405 


0.050 


0.410 


-0, 


.005 


35, 


,8 


326.4 


1.72 


54310.776 


0.477 


0.050 


0.454 


0, 


.023 


-31, 


.5 


326.8 


2.91 


54310.786 


0.558 


0.054 


0.619 


-0, 


.061 


-45, 


.5 


325.4 


3.16 


54310.797 


0.870 


0.100 


0.745 


0. 


.125 


-59, 


.5 


323.5 


3.42 



Table 3. Spectroscopic Orbital Solutions for CrB 



n ,1 t^TTl CiT\\ 


TViis; Wnrlr 


OIXVJO VCiUtjiLlCO 




Orbital Elements: 








P (days) 


1.139791423 ± O.OOOOOOOSO^^ 


1.139791423 ± 0.000000080'^ 


1.1397912 (adopted) 


T (HJD-2,400,000)'= 


50,127.61845 ±0.00020'' 


50,127.61845 ±0.00020^^ 


50,127.6248'' 


e 


O.O'^ 


O.O'' 


0.0<= 


w (deg) 


O.O'' 


0.0'' 


O.O'' 


7 (kms^"'") 


-13.03 ± 0.11 


-12.58 ± 0.05 


-12.3 ± 0.06 


ifp (kms-i) 


61.25 ± 0.21 


61.31 ± 0.06 


61.34 ± 0.06 


(kms^i) 


63.89 ± 0.22 


62.90 ± 0.08 


62.91 ± 0.08 


Derived quantities: 








MpSin^i [Mq) 


0.11818 ± 0.00092 


0.11461 ± 0.00032 


0.1147 


MsSin^i (Mq) 


0.11329 ± 0.00086 


0.11170 ± 0.00027 


0.1118 


q = Ms/Mp 


0.9586 ± 0.0047 


0.9746 ± 0.0016 


0.975 ± 0.002 


ftp sini (10^ km) 


0.9600 ± 0.0033 


0.96085 ± 0.00097 


0.96138 ± 0.00093 


flg sini (10® km) 


1.0014 ± 0.0035 


0.98592 ± 0.00126 


0.9861 ± 0.0012 


a sin I (Rq) 


2.8181 ± 0.0068 


2.7971 ± 0.0023 


2.798 ± 0.002 


Other quantities pertaining to the fit: 






Nobs 


46 


217 


217 


Time span (days) 


2610 


5.4 


5.4 


CTp (kms~-^)^ 


1.04 


0.74 


0.71 


(Ts (kms^"'^)^ 


1.10 


0.97 





^Our orbital solution using SR03 velocities. 
''Determined using all published velocities (see ij |2.ip 
'^T is the epoch of maximum primary velocity. 

''The value from SR03 has been shifted by an integer number of cycles to the epoch derived in this work. 



for comparison purposes. 
''Circular orbit adopted. 
fRMS residual from the fit. 



Table 4. Visual Orbit Solution for CrB 



Orbital Parameter Value 

Adopted values: 

Period (days) 1.139791423 ± 0.000000080=' 

Tnodc (HJD-2,400,000)'^ 50,127.04855 ± 0.00020 

e 0.0'= 

uj (deg) 0.0^= 

Op (mas) 0.50 ± 0.05^^ 

(mas) 0.45 ± 0.05^1 

Visual orbit parameters: 

a (mas) 1.225 ± 0.013 

i (deg) 28.08 ± 0.34 

n (deg) 207.93 ± 0.67<= 

/^K' 0.19 ± 0.19 

Reduced 0.61^ 



''See §0 

''This is the epoch of the ascending node, defined as the 
epoch of maximum secondary velocity, and accordingly is 
one-half period less than the value in Table [3] (see §[5]) . 

Circular orbit adopted. 

''See §E1 

''This value suffers from a 180'' ambiguity due to the co- 
sine term in Equation ([2]). 

^The low indicates that our error estimates for visibil- 
ity are conservative. 



Table 5. Physical parameters for cr^ CrB 





This Work 


SR03 Spectroscopy^ 


SR03 Results 


Physical Parameter 


Primary 


Secondary 


Primary Secondary 


Primary Secondary 


a (Rq) 


5.99 ± 0.07 


5.99 ± 0.07 




Mass {Mq) 


1.137 ± 0.037 


1.090 ± 0.036 


1.128 ± 0.037 1.099 ± 0.036 


1.108 ± 0.0041^ 1.080 ± 0.004*^ 


Radius (i?©) 


1.244 ± 0.050 


1.244 ± 0.050 


1.244 ± 0.050 1.244 ± 0.050 


1.14 ± 0.04 1.14 ± 0.04 


Teff (K) 


6050 ± 150 


5870 ± 150 


6000 ± 50 5900 ± 50 


6000 ± 50 5900 ± 50 


Mv (mag) 


4.35 ± 0.02 


4.74 ± 0.02 


4.45 ± 0.02 4.61 ± 0.02 


4.61 ± 0.07 4.76 ± 0.07 


Mk (mag) 


2.93 ± 0.09 


3.12 ± 0.11 







'^The se parameters use the SR03 spectroscopic results such as flux ratio, rotational velocities, and radial velocities, but 
use the iLestrade et al.l ()l999r ) parallax, Tycho-2 magnitudes, and our visual orbit. 

•^As noted in §[T1 these uncertainties are unrealistically small. 
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Table 6. Visual Orbit Solution for cr^-a'^ CrB 



Orbital Parameter 


Value 


P (years) 


726 ± 62 


To (BY) 


1,825.2 ± 1.5 


e 


0.72 ± 0.01 


w (deg) 


237.3 ± 6.8 


a (arcsec) 


5.26 ± 0.35 


i (deg) 


32.3 ± 4.1 


^ (deg) 


28.0 ± 0.5 
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Fig. 1. — Our radial velocities and the orbital fit for cr^ CrB (top panel) and the primary and 
secondary residuals (bottom panels). Filled circles represent the primary and open circles 
represent the secondary component. The corresponding orbital elements are listed in Table [31 
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Orbital Phase 

Fig. 2. — Same as Figured] but based on SR03 radial velocities. 
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Fig. 3. — Residuals for the individual nights' velocities from SR03. 
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Fig. 4. — Calibrated visibility measurements for cr^ CrB versus the projected baseline. The 
plus signs are the calibrated visibilities with vertical error bars, and the diamonds are the 
calculated visibilities for the best-fit orbit. Table [2] lists the numeric values corresponding to 
this plot. 
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Fig. 5. — The visual orbit of CrB. Open circles mark the positions of the two components 
at the epoch of ascending nodal passage, and the X marks identify the secondary's calculated 
positions at the epochs of visibility measurement. 
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Fig. 6. — distribution around the best-fit solution for the angular semimajor axis (a). The 
bottom dashed line corresponds to the minimum value, and the others mark a deviation 
of 1, 4, and 9 units above the minimum, corresponding to 1, 2, and 3-a errors. 
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7. — Same as Figure [6l but for the orbital inclination (i). 
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Fig. 8. — Same as Figure [6l but for the longitude of the ascending node 
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Fig. 9. — Same as Figure [6l but for the i^'-band magnitude difference {AK'). 
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Fig. 10. — The position of the Sun-hke components of a CrB on the H-R diagram. The points 
from top to bottom arc CrB primary, CrB, and o"^ CrB secondary. The isochrones are 
from the Yonsei-Yale (dotted) and Victoria-Regina (dashed models) for 0.5, 1.5, 3.0, and 5.0 
Gyr ages (left to right) for solar metallicity stars. 




Fig. 11. — The visual orbit of the wider cr^ — CrB (AB) system based on all measures 
in the WDS. Plus signs indicate micrometric observations, asterisks indicate photographic 
measures, open circles indicate eyepiece interferometry, and filled circles re present s peckle 
interferometry. The solid curve is our orbit fit and the dashed curve is the IScardial (Il979l ) 
orbit. — C lines connect each measure to its predicted position along the orbit. The big 
plus at the origin indicates the position of the primary and the dot-dash line through it is 
the line of nodes. Scales are in arcseconds, and the curved arrow at the lower right corner 
by the north and east direction indicators shows the direction of orbital motion. 
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Fig. 12. — Relative separation between CrB and ADS 9979C based on 10 resolutions of 
the pair from 1832 to 1984. Plus signs indicate micrometric observations. O—C lines connect 
each measure to its predicted position along the linear fit (thick solid hne) . The thick dashed 
line is the predicted movement based on the differential proper motions. The long dashed 
line connected to the origin indicates the predicted closest apparent position. The scale is in 
seconds of arc. An arrow in the lower right corner by the north and east direction indicators 
shows the direction of motion of the star. 
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Fig. 13. — Same as Figure fT2l but for ADS 9979D based on 106 resolutions of the pair from 
1825 to 1996. Asterisks indicate photographic measures and filled circles represent Tycho 
measures. 
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Fig. 14. — Mobile diagram of a CrB and some of its properties. The Ca-Cb pair is WDS 
component E, while WDS components C & D are not gravitationally bound to the a CrB 
system (see Figures [12] and [131 and text in §[7]). Op for the Ca-Cb pair is the photocentric 
semimajor axis. 



